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Abstract 

With use of CompHEP package we've made the detailed estimate of the influence 
of double e^e~ pair production (DPP) by photons on the propagation of ultra high 
energy electromagnetic (EM) cascade. We show that in the models in which cosmic 
ray photons energy reaches few x 10^ EeV refined DPP analysis may lead to substantial 
difference in predicted photon spectrum compared to previous rough estimates. 

1 Introduction 

Ultra-high energy (UHE) photons have not been recognized so far by any of present gen- 
eration experiments [H [21 El H] , although their existence is predicted by Greisen-Zatsepin- 
Kuzmin effect [SI E] as well as by most of hypothetical top-down models of UHE cosmic rays 
origin. There are several bounds on fraction and flux of ultra-high energy photons above 
10 — 100 EeV obtained by independent experiments [3, El E]- Photon limits are used to con- 
strain the parameters of top-down models (see for example [IH])- Future bounds may also 
limit considerable part of parameter space of astrophysical models, in which photons are pro- 
duced as secondaries from interactions of primary protons or nuclei with cosmic microwave 
background (CMB). Understanding interactions of UHE photons with universal backgrounds 
is a crucial point for building such constraints. 

In the wide energy range the spectra of electron and photon components of cosmic rays 
follow each other due to relatively rapid processes transferring 7-rays to electrons and back- 
wards. Pair production (PP) and inverse Compton scattering (ICS) are the main processes 
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that drive the EM cascade. In the Klein-Nishina hmit where s 3> ml, either electron or 
positron produced in a pair production event typically carries almost all of the initial total 
energy. The produced electron (positron) then undergoes ICS losing more than 90 % of 
energy and finally the background photon carries away almost all of the initial energy of 
the UHE photon. Due to this cycle the energy loss rate of the leading particle in the EM 
cascade is more than one order of magnitude less than interaction rate. However in presence 
of a random extragalactic magnetic field (EGMF) the electrons may lose substantial part of 
their energy by emitting synchrotron radiation. In this case, starting from certain energy 
the synchrotron loss rate for the electrons becomes to dominate over ICS rate, which leads 
to suppression of the EM cascade development. Its penetration depth is then defined by the 
photon mean free path. Depending on the value of EGMF this transition may occur between 
~ 1 EeV and ~ 10^ EeV. 

In this article we consider higher order process, double e~^e~ pair production (DPP) by 
photons. The DPP cross section grows rapidly with s near the threshold and quickly ap- 
proaches the asymptotic value a{oo) ~ 6.45 fih [HI [121 113] . The explicit energy dependence 
of the DPP cross section was estimated in Ref. [H] by calculating the dominant contribution 
from two e^e~ pairs to the absorptive part of gamma-gamma forward scattering amplitude. 
Since PP cross section decreases with the increase of ^/s the DPP rate starts to dominate 
over PP rate above certain energy. For interactions with CMB the transition occurs above 
~ 1000 EeV. In presence of the radio background this energy goes up somewhat. If the 
EGMF is less than 10^^^ G the EM cascade still exists at these energies and one should 
accurately count the secondary electrons from DPP. So far the EM cascade simulations such 
as [151 IB] roughly estimated DPP effect by utilizing the total cross section and assuming 
that one e~^e~ pair of the two carries all the initial energy while two particles in the pair 
are produced with the same energy. By making use of CompHEP package [T71 [THl [IS] we 
numerically calculate differential cross section for DPP and compare the influence of DPP 
on propagation of ultra high energy EM cascade with previous estimates. 

The paper is organized as follows. In Sec. [2] we present the results of the calculation of 
DPP cross section. In Sec. |3] we write transport equations for EM cascade and calculate 
the coefficients for transport equations for photons, and secondary e~^,e~ related to DPP. 
In Sec. IHwe illustrate the influence of DPP in model example. In Sec. |5]we summarize our 
results. 
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2 DPP cross section calculation 



As it was mentioned in the Introduction the DPP process begins to dominate over PP at 
very high energies > 1000 EeV or s > 1 GeV^ which is well beyond the DPP threshold. 
At these and higher energies DPP has noticeable effect on the propagation of EM cascade. 
In this energy region DPP total cross section is practically saturated by its asymptotic value. 
So, we are interested here mostly in energy and angular distributions of secondary electrons 
(positrons) in asymptotic regime {s oo). 

We use CompHEP package for calculation of tree level differential DPP cross sections. 
This package allows to perform automatic calculations of matrix elements and their squares 
for any process 2 — 2, .., 2 ^ 4 at tree level. Then, with the aid of CompHEP one 
can integrate squared matrix elements over selected part of multi-particle phase space. See 
Refs. HH IIHI [IS] for the details. 

We introduce binning in the energy E* of one of the produced electronic]. Then we 
perform CompHEP simulations in the centre of mass frame (CMF) and obtain distributions 
over cosO* of the cross section in a given energy bin. Here B* is the angle between the 
collision axis and the momentum vector of the electron. 

These calculations show that the angular distribution of secondaries tends to a strongly 
peaked function of cos G* in the asymptotic energy range. The peaks are located at forward 
and backward directions, i.e. at cos0* = ±1. This behavior is illustrated by Fig. [T] where 
we show an example of angular distribution for y/s = 10.0 GeV. The effect that most of 
secondaries go forward or backward becomes more pronounced with the increase of ^/s and 
in the case of fixed ^/s with higher energy E* of electron. Numerically we found that the 
probability of emitting secondary electron inside the cone with | 1 — cos©* |< 1/50 is 96.8% 
for = 1.0 GeV, 98.7% for = 2.5 GeV and 99.6% for = 10.0 GeV. Also we checked 
that the probability of producing two forward secondaries of the same type (e.g., when both 
forward particles are electrons) integrated over energies and directions of other secondaries, 
is of order 10~^. So, the main part of events consists of two e~^e~ pairs going to the opposite 
directions along the collision axis. 

Let us now turn to the energy distribution. It is clear from the symmetry of the problem 
that the energy distribution in the CMF frame should be the same for the forward and 
backward electrons. Let us write the DPP differential cross section in the form 

^Here and further we denote by "*" quantities measured in the center of mass frame. 
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Figure 1: Example of distribution of the DPP cross sections over cos0* for ^/s = 10.0 GeV 
in energy E* bin 2.25 — 2.5 GeV. 



Then the energy conservation condition gives 

da 

VscTtotis) =4: J ^*'^ = ^tot{s) Vs J r(f){r, s) dr 
or for any value of s 

I r(j){r^ s)dr = 1 . (2) 
Jo 

Imposing probability conservation requirement 



da 







gives another integral constraint on 0(r, s): 

I (f){r^ s) dr = 2 . 
Jo 



(3) 



Although conditions ^ and ^ do not necessary imply d(f)/ds = the results of CompHEP 
simulations show that for large enough s, when cross section approaches its asymptotic value, 
the energy distribution of secondaries in units of maximal energy v^/2 varies only slightly 
with a/s. In Fig.[2]we plot the distribution 0(r, s) as a function of r for different values of a/s. 
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Figure 2: The normalized energy distributions (j){r, s),r = 2E* / ^/s (see main text for details): 
1) ^ = 0.1 GeV, 2) y/s = 0.25 GeV, 3) ^/s = 10.0 GeV and 4) the analytic approximation 
Eq. ©. 

One can see that with varying ^/s the only changes in these distributions are concentrated 
at the borders of the plot. The distribution limit for s ^ oo can be fitted (see fig. [2]) by 
simple analytic expression, which satisfies constraints and ([3]), 



For comparison, the earlier approximation [151 KB] in terms of the distribution 0(r, s) reads 



In our further calculations we use for the energy distribution the Eq. ([T]), where 0(r, s) = 
(pfitij) is given by Eq. (jlj) and assume that all the secondary particles are directed alongside 
the collision axis. In Sec. H] we discuss how good the above approximation is. 
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0(r,s) = 2-(5(r-0.5). 



(5) 



3 Transport Equations 



Here we describe propagation of the UHE cosmic rays using the formalism of transport 
equations in one dimension. Besides DPP term on which we are going to focus now, the 
full transport equations for the electrons and photons contain the terms describing ICS, PP, 



synchrotron and e~^e~ pair production by electrons and positrons as well as redsliift terms. 
For simplicity, we show below the part of the equation written for nonexpanding universe 
with the terms related to the DPP process only. 

-N,{E,,t)= dE,N,{E„t) den{e) dfi^^^{E,,E^,s) (6) 

j^N^{E^,t) = -N^{E^,t) f den{e) I dfi^^a^ppis) (7) 



where Ne{Ee,t) is the (differential) number density of electrons at energy E^ at time t, n{e) 
is the number density of background photons at energy e, fi is the cosine of the collision angle 
(/i = — 1 for a head-on collision) and s = 2E^e{l — fi) is center of mass energy squared. The 
term in the r.h.s. part of Eq. ([U]) describes influx of electrons produced in DPP. Transport 
equation for positrons has the same form as The r.h.s. term of ([7]) describes the loss 
of photons due to DPP. The factor (1 — /i)/2 is the flux factor. As we've seen in the Sec. [2] 
the pairs produced in DPP are directed alongside the collision axis. This implies that one 
of the pairs carries practically all the initial energy of the photon in the laboratory frame. 
Here we neglect the nonleading pair produced in the interaction. 
Replacing integration over fi by integration over s gives 



where 



h{x) = / -^de . (10) 



Here Sth = IQml is threshold CMF energy squared for DPP and Smax = ^E^emax- 

Now we are ready to use the results obtained in the Sec. [21 Again here we calculate 
the transport equation coefficients in the limit of s^Sth- This implies that electrons and 
positrons are ultrarelativistic in the CMF frame. The CMF 7-factor in the laboratory frame 
is 

1 E 

ICMF = {l-PhMp)'^ = -r ■ 

Provided that e"*" and e~ momenta are directed either towards the CMF frame velocity or 
in the opposite direction, their energy in the laboratory frame 

Ee = 7cMFE:ii ± /?:) = ^E:ii ± /?:), 
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where /3* — >■ 1 is electron velocity in CMF. For the leading e^e pair we have: 
Then using Eq. ([T]) we finally obtain 

Using numerical simulations of cosmic rays propagation presented in the Sec. lUwe have also 
verified that utilizing simple step function for the total cross section 

o-DPp(s) = crDpp(oo)e(s - Sth) (12) 

instead of exact one listed in ^14j, doesn't introduce any visible change to the resulting 
spectra. This implies that the equations ( fTTj) and ([9]) can be simplified as follows: 

^N4E,,t) ^ aDPp(oo) r dE,^^&^(j)&K,{sth/^E,), (13) 
at 

j^N^iEy, t) ~ -2(TDPp(oo)Ar^(E^, t)K,{sth/^E^). (14) 

where 

K,ix)= / iMydy (15) 

J X 

is the function totally determined by the background photons spectrum. 



4 Model example 

In the previous sections we have found the precise expression for the distribution of secondary 
electrons from DPP. Here we consider a model example to illustrate the difference introduced 
by the specified cross section compared to the previous estimates. 

We use a numerical code developed in Ref. [13] to compute the fiux of produced photons 
and protons. The code is based on the transport equations and calculates the propagation of 
nucleons, electrons and photons using the dominant processes. For EM cascade it includes 
all the processes mentioned above. For nucleons, it takes into account single and multiple 
pion production and e~^e~ pair production, neutron /3-decay. The propagation of nucleons 
and the EM cascades are calculated self-consistently, that is secondary particles produced in 
all reactions are propagated alongside the primaries. 

Besides CMB the radio, infra-red and optical (IRO) components of the universal photon 
background are taken into account in the simulation. Note that the radio background is 
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not yet well known. Our results will depend strongly on the radio background assumed. 
Three models considered in this work are estimates by Clark et al. [20] and the two models 
of Protheroe and Biermann [21j, both predicting larger background than the first one. For 
the IRO background component we used the model [22j. This component doesn't have 
substantial effect on the propagation of UHE protons and EM cascade. For the strength of 
the random extragalactic magnetic field we use the range of values lO^^^G < B < lO^^^G 
following the estimate [23] . 

Among the models we have chosen the one in which the UHE photons contribute sub- 
stantial part of the total spectrum. Note that such models are strictly limited by the present 
experimental bounds on the photon component, see [TU] for details. 




0.5 ■ 
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Figure 3: Fluxes predicted by proton emitting source described in text. The dotted line 
represents proton component, while solid, dashed and dash-dotted lines represent total fiux 
calculated with CompHEP based DPP, trivial DPP estimate and without DPP correspond- 
ingly. 



In Fig. [3] the propagated cosmic ray fiux is shown for proton sources with spectrum 

^ ~ ■^ E < lO^EeV (16) 
dE 

homogeneously distributed in the Universe and having no evolution in the comoving frame. 
The spectrum presented is normalized on HiRes [3] results (fitting was done above 40 EeV). 
The solid line represents the total UHE cosmic ray fiux calculated with use of new DPP esti- 
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mate. The dotted line shows proton component. The dashed hne shows total flux calculated 
using earlier DPP estimate (Eq. [3]), utilizing the total cross section and assuming that one 
e^e~ pair of the two carries all the initial energy while two particles in the pair are produced 
with the same energy. The dash-dotted lines are built without taking DPP into account at 
all. 

It is clear from the Fig. that DPP suppresses 7 ray flux above 100 EeV. This is only 
true if the minimal radio background model |20j is used. The same picture made for any 
of the two models of [2Tj haven't shown any effect of DPP, since in this case the 7 flux is 
strongly suppressed by PP on radio. Increasing magnetic field above lO^^^G also destroys 
the picture, this time due to synchrotron radiation. In the case of minimal radio background 
and moderate EGMF the trivial DPP effect estimate leads to extra suppression compared 
to the more accurate one proposed in this paper. Although overall error in terms of integral 
photon flux above 100 EeV turns to be only +7% for the curve disregarding DPP and just 
— 1.5% for the trivial DPP estimate. Note that integral photon flux fraction predicted in 
this model is 34%, which is very close to the upper bound [7]. So far we used the fixed 
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Figure 4: Electron flux predicted by photon emitting source described in text using 0) earlier 
estimate Eq. ([5]); 1) analytical fit (jlj); 2) analytical fit (jl]) + 3.2% perpendicular component 
(see details in text); 3) 0(r, 1 GeV^); 4) 0(r, 100 GeV^). 

energy distribution (jll) for all values of s. Also we assumed that all the secondary particles 
are directed alongside the collision axis. Let us now check how accurate the above approxi- 
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mations are. We have repeated our simulations replacing the energy distribution (jl]) by the 
tabulated functions obtained with use of CompHEP for ^/s = 1 GeV and = 10 GeV. To 
see the maximal possible effect of the nontrivial angular distribution of secondaries we have 
also repeated our calculations assuming that 3.2% of secondary particles are aligned perpen- 
dicular to the collision axis in the CMF, while the rest of the particles are directed alongside 
the axis. We don't show here the modified fluxes obtained in the model corresponding to 
Fig. [3l since they are practically indistinguishable from the curves already shown. Instead 
to illustrate the maximal possible error introduced by the approximation used, here we con- 
sider the pure photon sources with the same injection spectrum (I16p as in Fig. [H] and count 
the income to the propagated electron and photon spectra from the uniformly distributed 
sources located within 500 Mpc from the observer. In Figs. H] and [5] the electron and photon 
fluxes in this model are shown respectively. Also on these figures the fluxes calculated using 
10 1- 

9 - 

8 - 
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Figure 5: Propagated photons flux predicted by photon emitting source described in text. 
The designation of curves is the same as on Fig. |H 




earlier estimate Eq. ([5]) are shown. From the figures it is clear that the earlier estimate may 
lead to the artificial features in the spectra which doesn't appear in our analysis. Also it is 
clear that discrepancy between the curves 1-4 representing different variants of our analysis 
are small compared to the error introduced by the earlier estimate. In fact the difference 
between the curves 1-4 is comparable to the error introduced by finite energy binning used 
in our numerical code. 
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5 Conclusion 

In this work we have considered in detail the DPP process. We have estimated the distri- 
bution of secondary electrons and positrons and made the improved cosmic rays simulation 
based on the new estimate. We have shown that in certain cases the DPP process may 
modify the photon component of the spectrum substantially. However this modification can 
only be seen if radio background is close to the minimal model [20] and EGMF is lower than 
lO^^^G. In this case there is an energy range where DPP is the main attenuation mechanism 
for 7 rays and therefore differences in DPP estimates can clearly be seen. Although in the 
vast majority of the models which do not contradict to the present experimental bounds on 
the photon fraction in UHE cosmic rays DPP process doesn't make a substantial contribution 
to the attenuation and therefore can be treated simplistically or even disregarded. 
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